Penetration depth anisotropy in two-band superconductors 
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The anisotropy of the London penetration depth is evaluated for two-band superconductors with 
arbitrary inter- and intra-band scattering times. If one of the bands is clean and the other is dirty in 
the absence of inter-band scattering, the anisotropy is dominated by the Fermi surface of the clean 
band and is weakly temperature dependent. The inter-band scattering also suppress the temperature 
dependence of the anisotropy. 
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The two-gap superconductivity of MgB2 is estab- 
lished experimentalljiiiSiMi^ and by solving the Eliash- 
berg equations for the gap distribution on the Fermi 
surfaced According to the latter, the gap on the four 
Fermi surface sheets of this material has two sharp max- 
ima: Ai w 1.7 meV at the two 7r-bands and A 2 ~ 7meV 
at the two cr-bands. Within each of these groups, the 
spread of the gap values is small. 

A number of physical properties of MgB2 were rea- 
sonably well described within a model with two con- 
stant gaps on two separate Fermi sheets. Still, the 
data on anisotropy of the magnetic field penetration 
depth A are controversial. The anisotropy parameter 
7a = A c /A a has been calculated within the weak- coupling 
clean-limit model and shown to increase from about 1 . 1 at 
T = to ss 2.6 at T c & Similar prediction has been made 
within Eliashbcrg formalism. 9 Qualitatively, the predic- 
tions were confirmed in STM^fiiAi small angle neutron 
scattering^ and magnetization experiments pi£ However, 
other groups recorded different behavior Jitifiiik Given 
variety of samples used, it seems imperative to consider 
effects of scattering upon 7a, a non-trivial problem given 
different roles of the intra- and inter-band scattering in 
two-band materials. 

In the following, we reiterate that in the presence of 
inter-band scattering, the energy gap in the quasiparticle 
excitation spectrum as revealed by the density of states 
(DOS) differs from both A1 and A 2 r^ the situation 
reminiscent of the Abrikosov-Gor'kov pair breaking. 19 
For this reason, we use the term "order parameter" for 
A's. We stress that all thermodynamic properties depend 
on the actual DOS and are affected by the inter-band 
scattering. Then, we show that 7a = A c /A a {, depends on 
both inter- and intra-band scattering. 

Our approach is based on the quasiclassical version of 
the BCS theory for a general anisotropic Fermi surface 
and for an arbitrary anisotropic order parameter A(fc). 20 
In the absence of currents and fields we have for the Eilen- 
berger Green's functions f(k,u>) and g(k,u>): 

= 2A 5 -2a;/ + 7, 1 = g 2 + f 2 . (1) 

Here the scattering term I is given by the integral over 
the full Fermi surface: 

I(k) = [ d 2 q p(q) W(k, q) [g(k)f(q) - f(k)g(q)} (2) 



with W(k, q) being the scattering probability from q to 
k at the Fermi surface. The Matsubara frequencies are 
lu = nT{1n + 1) with an integer n {% = 1). The local 
DOS p{q) is normalized: J d 2 qp(q) = 1. 

This system 0-10 should be complemented with an 
equation for A(fe). We will not use it here, rather tak- 
ing A(fe) as a given. This simplifies the problem greatly 
because solving for A(fe) usually involves a number of 
assumptions which are difficult to control. 

We use the approximation of the scattering time r: 

Jd 2 qp(q)W(k,q)f(q) = (f)/r; (3) 

(...) stands for the average over the Fermi surface. 
Clearly, the approximation amounts to the scattering 
probability W = 1/r being constant for any k and q. 

For two well-separated Fermi surface sheets, the prob- 
abilities of intra-band scatterings may differ from each 
other and from processes involving k and q from different 
bands. The effects of the inter- and intra-band scattering 
upon various properties of the system are different, e.g., 
the intra-band scattering does not affect T c , whereas the 
inter-band does. Therefore, Eq. is replaced by 

J d 2 q a p(q a ) W{k/3, q a ) f{q a ) = (f) a /rp a . (4) 

Here a,(3 = 1, 2 are band indices and (...) Q denotes aver- 
aging only over the a-band. 

We now assume the order parameters A(fe a ) taking 
constant values Ai and A 2 at each of the two bands. 
Writing Eq. for k in the first band, we have: 

=2A 1 g 1 -2uf 1 

+ — [9i{fh-h{9h} + — [<?i</>2-/i<<?> 2 ]. (5) 

Til Ti2 

For a uniform sample in zero field and with k independent 
A's in each band, the functions f,g are k independent: 
(f)a = fa and (g) a = g a . Then, we have: 

= Ai#i - wft + [51/2 - /i02]/2tl2 ■ (6) 

The equation for the second band differs from this by 
replacement 1 «-> 2. The fact that rn and T22 do not 



2 



enter the system (JHJ is similar to the case of one-band 
isotropic material for which non-magnetic scattering has 
no effect upon T c (the Anderson theorem). It is the inter- 
band scattering that makes the difference in the two-gap 
case, the fact stressed already in early workiiZiiS 

Two equations JSJ) are complemented with normaliza- 
tions + f 2 = 1 to form a sufficient set. Following Rcf. 
Il9l we introduce variables u a = g a /f a and obtain after 
simple algebra*iii2i2i 



a7 = Ui + Ci 



u 2 



A 2 



— = U 2 + Ci 



u 2 - Ul 



Ci = 

(2 = 



1 



2t 12 A 1 
1 

2r 2 iA 2 



The Eilenberger functions in terms of variables u read 
/ = 1 



(7) 



1 



g = u 



1 



(8) 



In general, the system (jjTj) can be solved only numeri- 
cally. However, near T c , u — g / f 1 and one obtains: 



Ml = 



uj uj + CiAi + C2A2 
A7 oj+ (Ci +C2)A 2 



(9) 



w 2 is obtained by 1 <-* 2. Clearly, u Q = uj/A. a in the 
absence of inter-band scattering. For £ ^> 1, we have: 



Ul « W 2 



e* 



(Ci +C2)AxA 2 
CiAi +C 2 A 2 



(10) 



Moreover, if the inter-band scattering is strong, Eq. 
(|10|l holds at any T. To see this, look for solutions of 
Eqs. Q) in the form 



UJ 



1,2 



(11) 



where v a are small corrections. Substitute these in Eqs. 
@ and keep only linear terms in v to obtain 



ff*(e*-A!) 
A 1 [l + 5 *(Ci+C 2 )] 
e*(CiAi + C 2 A 2 ) - AiA 2 (Ci + C2) 
A 1 A 2 [l+ 5 *(Ci+C 2 )] 



(12) 



where g* = uj/\/uj 2 + e* 2 . For £ Q — > 00, v\ remains small 
only if e* is given by expression (|10|l . 

The DOS for two-band materials is N(e) = 
N(0) Re(^i gi + v 2 g 2 )u>^ie] v a are fractions of the to- 
tal DOS N(0) at the two pieces of the Fermi surface. 
For strong inter-band scattering, this gives in the lowest 
approximation 



N(e) = N(0) 



VP 



(13) 



i.e., e* is the common for both bands energy gap. 

It does not seem possible to provide a general expres- 
sion for e* in terms of A Q and an arbitrary inter-band 



scattering strength. Still, in principle, we can evalu- 
ate any thermodynamic property of a two-band material 
knowing the solutions u of the system (J2J and the gap e* . 

If the ground state functions (which we call now /(°), 
<?(°)) are known, one can study perturbations of the uni- 
form state such as penetration of a weak magnetic field, 
i.e., the problem of the London penetration depth. The 
perturbations, /W, g( l \ should be found from the full 
Eilenberger equations^ we have for the first band: 

vTlfi = 2A m - 2uh + — [,gi(/)i - A (.9)1] 

Til 

+— M/) 2 -/i<s) 2 ], (14) 

T~12 

Here, v is the Fermi velocity, II = V + 2iriA/<fio. The 
second equation is obtained by 1 <-> 2. Two equations for 
the "anomalous" functions / + are obtained from these by 
complex conjugation and by v — ► — vM' The normaliza- 
tions g 2 + f a f£ — 1 complete the system. 
We look for solutions in the form 



fa = (fL 0) + fL 1} )e^ r \ ft = (fL 0) + fP + )e- i6 ^ , 
9c = 9^+9^, a = 1,2. (15) 

where f a °^ and g a °^ can be expressed in terms of u's ob- 
tained solving the system Q). The form (f 1 51) takes into 
account that in the London approximation only the over- 
all phase depends on coordinates. We obtain for the cor- 
rections after some algebra: 



ff «Ai-/WV = iA (0, «i'/2, 



p(0) 



(16) 



2 5 iV ) +/rar+/r") = o, 



e(0){Al) 



where P = V# + 2ttA/(/)q and 

A; = Ax + / 1 (0) /2r 11 + / 2 (0) /2n 2 , (17) 

wj = uj + gT /2r n + g ( 2 0) /2r 12 . (18) 

The equations for the second band (decoupled from the 
first) are obtained by 1 «-> 2. 

To evaluate the penetration depth we turn to the Eilen- 
berger expression for the current density^ 



j = -47r|e|iV(0)T Im^(vg) 
and compare it with the London relation 

47T 



-^-(A 2 )^(|> + A 



(19) 



(20) 



Here, (A 2 )^. 1 is the tensor of the inverse squared pene- 
tration depth (proportional to the superfluid density ten- 
sor); summation over k is implied. We now find g± from 
the system (tTBll : 



9? 



■ A0)2 n 
lf{ VP 



2(A' 1 / 1 (0) +ojW< 



(0)- 



f (0)2 (0) 

= i^vP; (21) 
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<72 i s obtained by replacement 1 <-> 2. Substituting g a ' 
in Eq. (|19fl and comparing with Eq. (|20[) we obtain: 

16^ 2 e 2 7V(0)T ^ , . f 2 g a 



(i) 



(22) 



at two limiting situations when one of the bands is clean 
while the other is dirty. If the first (tt) band is clean and 
the second (cr) is a dirty extreme (T22A2 — ► 0), one can 
disregard the contribution of the dirty band to obtain for 
both T = and T = T r : 



Only the unperturbed functions /, g enter the penetra- 
tion depth; for brevity we dropped the superscript (0). 
Equation \\l'2\i is our main result. Thus, to evaluate the 
penetration depth for given order parameters A a in the 
presence of scattering, one has to solve the system (J7J) 
for u a (u>), then to substitute the equilibrium functions 
fa, ga (given in ©) in Eq. to sum up over uj. 

The band calculations^ yield for MgB2 the following 
averages: (v 2 )i = 33.2, (v 2 )i = 42.2, {v 2 a ) 2 = 23, and 
(v 2 }2 — 0.5 x 10 14 cm 2 /s 2 . Tensors (viVk)i and (v l Vk)2 
have opposite anisotropics: 



(23) 



whereas averaging over the whole Fermi surface yields a 
nearly isotropic result: (v 2 )/(v 2 ) w 1.2 . 

In the clean limit (all r a p — > 00) u>' — uj and A' a = A a . 
Besides, u a = uj/A a and f 2 g a /u' a = A 2 /(w 2 + A 2 f/ 2 . 
Expression i|22|) reduces to the result given in Rcf. 
0. For MgB2, it gives nearly isotropic penetra- 
tion depth at low temperatures: at T = the 
sums over uj in Eq. (|22|l are the same; this gives 

7a(0) = Xcc/Ka = « LI- Near T c , the 

sums are oc A 2 , and the contribution of the strongly 
anisotropic cr-band with the large gap dominates; this 
gives 7a (Tc) ~ 2.6. The curve 1 in Fig. ^ shows 7a (T) 
for this case. 



7a(0)= 7 a(T c ) 




0.89. 



(25) 



If the (tt) band is dirty and the (cr) is clean, we have 



7a(0) = lx(T c ) 




(26) 



These two numbers constitute the minimum and max- 
imum possible values for A-anisotropy of MgB2. Thus, 
when one of the bands is clean and the other is dirty we 
expect a weakly T dependent 7a, the value of which is 
determined by the clean band. 

If the intra-band scattering is strong in both bands 
(tiiAi ~ T22A2 <C 1, T12 = 00), the bands contribute to 
the superfluid density tensor (A 2 )^ 1 as two independent 
dirty superconductors. To see this, we note that lu[ « 
gi/2rn and the sums over to in Eq, (|22|l can be evaluated 
exactly: 



\ - fhh 



E 



2r n A 2 
w 2 + A? 



T11 Ai . , Ai 

tann — 

2T 2T 



(27) 



Then, we arrive at the result obtained by Gurevich with 
the help of the dirty limit Usadel equations* 2 ^ 



(^^E^^^ ( 28 ) 



c 2 h 



2T 



Zero inter-band scattering (t\i — ti\ —00). If only the 
intra-band scattering is present, the functions /, g are the 
same as in the clean limit. We readily obtain: 



flg a 



A 



P 2 a (P a + l/2r aa ) 



(24) 



with /3 2 



+ A 2 . This expression appears in the 



standard penetration depth calculations, see e.g. Rcf. 
I23L For known A a (T), the sums in Eq. J52J) can be 
evaluated numerically; however, for T — > 0, T c , and in 
the dirty limit they can be done analytically. 

At T — 0, the sums are replaced with integrals 
according to 2ttT^2 u = J °° dm. Denoting I(T) = 
27rT]T w A 2 /(3 2 ((3 + 1/2t) we obtain 1(0) = 1 for r -> 00 
and 1(0) = tttA for rA < 1. 

Near T c , g— >1 and we have for clean bands 
^ f 2 /io' = 7C(3)A 2 /8tt 3 T c 3 , whereas for dirty bands it 
is tA 2 /4T 2 . 

Different impurities introduced to MgB2 may affect 
differently the scattering within the bandsi 2 ^* 2 ^ It is of 
interest to see how the anisotropy of A is affected by dif- 
ferences in scattering times th and T22. We first look 



where the anisotropic conductivities of the two bands 
ajfc — 2e 2 (viVk) a TaaV a N(0) are introduced (we write 
here h explicitly to avoid confusion in dimensions). This 
yields: 



t!(o) 



lliTc) 



Cri 1 Q ) A 1 +Cri 2 a ) A 2 

a ( cJ A 1 + A 2 



crgAf + 4gA| 
cr^Af + cr^Al 



(29) 



(30) 



Finally, we discuss the possibility of a strong inter-band 
scattering. As was shown above, in this case u — ui/e* + 
O(l/0 in both bands, see Eqs. JTTJ and GD- Then, the 
Eilcnbcrgcr functions are also the same in the two bands: 

Evaluation of the 



/ = e*/v^ 



+ e* 2 , g = uj/\/uj 2 + ■ 



-*2 



sums over to in Eq. (|22|l is now simple: 



/a2\-i 87r 2 e 2 iV(0) _ , e' ^ 
( A )ik = 6 tanh 2T E u cx{viV k ) a T a (31) 
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If all t's are the same, we have: 



0.15 



0.10 



0.05 



o.oo 




T/T 



FIG. 1: The anisotropy 7a = A c A a6 and the inverse 
square of the penetration depth L 2 /\\ b versus T/T c ; L 2 — 
16n 2 e 2 N{0)(vl)/c 2 . The curves labelled 1 is the clean limit, 
all 1/t are zero. The curves labelled 2 and 3 are calculated 
for a weak inter-band scattering: ri2Ai(0) = 500, T21 A2(0) = 
2000 (h = 1); the curve 2 is for a clean 7r-band, rnAi(0) = 10, 
and a dirty a, T22A2(0) = 0.1; the curve 3 is for a dirty ir 
and clean a: rnAi(0) = 0.1, T22A2(0) = 10. Curves 4 are for 
the intermediate inter-band scattering strength Ti2Ai(0) = 5, 
t 2 iA 2 (0) = 20, and mAi(0) = 0.05, r 2 iA 2 (0) = 2. 



where e* (T) is given in Eq. (fTUf> and 



n 



T22T21 



(32) 



Thus, all components of (X 2 )ik have the same T depen- 
dence and the anisotropy parameter is T independent: 



7a 



^l(^q)l T l +^2(^a)27-2 
Vl{vf)m + V2(v 2 ^)2T2 



(33) 



ll = (v 2 a)/(v 2 c ) ■ (34) 



For T — ► T c , this result was obtained in Ref. |8| we now 
see that it holds at any temperature. 

To recover the behavior of 7a (T) between and T c one 
needs explicit dependencies A(T). Qualitatively, this be- 
havior can be studied assuming A Q (T) A Q (0)(1 — t 2 ) 
with, e.g., A 2 (0) = 4Ai(0) = 2T C . Figure 1 shows results 
of numerical evaluation of 7a (T) for scattering parame- 
ters given in the caption (which are not that extreme as 
in the above discussion) . The curves 7a (T) are obtained 
by solving Eqs. Q for it's in two bands and then by eval- 
uation of the sums in Eq. lll'L'1) . It is worth noting that 
although the T dependences shown in the figure are ob- 
tained using approximate A(T), the end points of these 
curves at T = and T = T c are exact. 

We conclude that both the inter- and intra-band scat- 
tering affect strongly the superconducting anisotropy of 
two-band superconductors in general and of MgB2, in 
particular. If one of the MgB2 bands is dirty, the 
anisotropy is dominated by a cleaner band: 7a (T) is close 
to unity (and might be even less than 1) if the n band 
is in the clean limit, whereas in the opposite situation of 
a clean a, 7a (T) is large being in both cases weakly T 
dependent. The inter-band scattering suppress the T de- 
pendence of 7a as compared to the clean limit discussed 
earlier A 
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